The subdiffusive target problem: Survival probability 
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The asymptotic survival probability of a spherical target in the presence of a single subdiffusive 
trap or surrounded by a sea of subdiffusive traps in a continuous Euclidean medium is calculated. 
In one and two dimensions the survival probability of the target in the presence of a single trap 
decays to zero as a power law and as a power law with logarithmic correction, respectively. The 
target is thus reached with certainty, but it takes the trap an infinite time on average to do so. In 
dimensions higher than two a single trap may never reach the target and so the survival probability 
is finite and, in fact, does not depend on whether the traps move diffusively or subdiffusively. When 
the target is surrounded by a sea of traps, on the other hand, its survival probability decays as a 
stretched exponential in all dimensions (with a logarithmic correction in the exponent for d = 2). 
A trap will therefore reach the target with certainty, and will do so in a finite time. These results 
may be directly related to enzyme binding kinetics on DNA in the crowded cellular environment. 

PACS numbers: 02.50.Ey, 05.10.Gg 



I. INTRODUCTION 

A number of recent experiments have pointed to the 
occurrence of subdiffusive motion in biophysical environ- 
ments [1 H H H i, i, 0, i, 11, that is, motion where 
the mean square displacement of a "walker" grows sub- 
linearly in time, e.g., 



(r 2 (t))~t~<, 0< 7 <1. 



(1) 



In particular, recent experiments on transport of large 
molecules in living cells indicate that the crowded cellu- 
lar environment leads to such motion [l|, |2j, y, |J, l5[ ■ The 
sublinear growth of the mean square displacement comes 
about because the presence of a large number of macro- 
molecules leads to a medium that may present impedi- 
ments such as barriers, traps, or otherwise interrupted 
pathways to the diffusive motions observed in typically 
much sparser environments designed in vitro. A num- 
ber of experimental and theoretical papers have begun 
to address the problem of how to characterize the mo- 
tion of proteins in living cells and, having characterized 
it as subdiffusive, of how to estimate the effects of this 
slowed-down transport on the binding and unbinding of 
enzyme proteins on DNA target sites d, [T3, H3 ■ The 
results are not obvious and involve the balance between 
reaching a target site more slowly (which might tend to 
slow down reaction rates) while at the same time retain- 
ing the enzyme in the vicinity of a target site for a longer 
time (which might tend to speed it up). Furthermore, 
the breaking of ergodicity indicates that one might need 
to be careful about dealing with time histories of sin- 
gle events vs ensemble averages over many events [Tl| . 
These issues arise not only in the biophysical context 
of cells but in other situations involving crowded envi- 
ronments such as porous media, but it is the cellular 
context that has awakened interest in the problem most 



recently. Various theories appear not to be entirely con- 
sistent with one another in their predictions. For exam- 
ple, Golding and Cox Q predict singular features in the 
binding/unbinding kinetics when the exponent 7 crosses 
the value 2/3, while others [l(J[ll[ do not seem to predict 
such features. 

The theoretical approaches to the problem are fairly 
complex and range from continuous time random walk 
formalisms to scaling arguments. The literature that fo- 
cuses on the biophysical context includes a number of 
realistic features such as the probability of a protein un- 
binding from the target site on the DNA before a re- 
action takes place and having to return to that site. 
These features add complications surrounding the appro- 
priate boundary conditions. The relation between differ- 
ent boundary conditions used in this problem [l(| E| are 
not entirely clear and seem not to have been discussed 
in the literature. Furthermore, it may be the case that 
the final results of interest obey scaling laws that relate 
the more realistic problem to the simpler one in which 
binding definitely leads to reaction In any case, it 
seems to us that even the "stripped-down" version of the 
problem wherein a reaction occurs with certainty when 
the target is reached would benefit from a systematic and 
clear presentation, and this is our goal in this paper. 

Given the scenario of a target surrounded by a sea of 
subdiffusive particles the first of which to reach the tar- 
get defines the binding time of interest, what exactly is 
it that one wishes to calculate? In some instances, the 
quantity calculated is the probability that a subdiffusive 
molecule initially at a distance r from the target ever 
finds it [f| [HI O EH ■ In other instances, it is the distri- 
bution of the time-averaged probability of the subdiffu- 
sive molecule to be in the bound state (vs the unbound 
state) in a single trajectory [TTj . A standard classic mea- 
sure is the distribution of first reaction times [l(| HH EH • 
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We focus on a traditional quantity from which many 
of these can in principle be deduced under appropriate 
physical conditions, namely, the surviv al p robability Q(t) 
of the target as a function of time [HI, [la, . A specific 
scenario for this calculation arises when a target (per- 
haps DNA) site is surrounded by many randomly located 
particles (perhaps site-specific DNA-binding proteins). 
From the outset, we focus on the problem in a spatial 
continuum and thus rely on the fractional diffusion equa- 
tion. We note that essentially all the dynamical theories 
on this topic in the literature either start from or arrive 
at a continuum formulation, even those that begin with a 
continuous time random walk, and we furthermore note 
that there seems to be little argument about the fact 
that at least asymptotically and on a spatial mesoscale 
such a formulation is appropriate. The differences in the 
literature lie not in the use of a fractional diffusion equa- 
tion but rather in the boundary conditions [l(J El, El • 
In this context we note that while a portion of the lit- 
erature deals with partially absorbing targets, that is, 
targets with a finite probability of letting the attached 
particle go (to eventually return and to again be trapped 
or not), or equivalently, with a finite rather than an in- 
finite reaction rate when a particle meets the target, we 
deal with the fully absorbing target. The connection be- 
tween the two types of targets may in many cases be 
straightforward [id ]. 

The paper is organized as follows. In Sec.[TT]we present 
the formalism for the calculation of the survival proba- 
bility. The results for all dimensions are presented in 
Sec. IIIII Finally, we conclude with a brief summary in 
Sec EH 



II. SURVIVAL PROBABILITY: THE 
FORMALISM 

We start by defining Qi(r, t; R) as the probability that 
a random walker ( "the particle" ) that started at location 
r at time t = has survived until time t in the presence 
of an absorbing sphere of radius R centered at the origin. 
Here r = |r| in recognition of the orientational symmetry 
of the problem. The probability Qi is the main quan- 
tity from which other results are obtained. If instead 
of one we have N independent (sub)diffusive particles, 
their combined survival probability Q{t\ R) is simply the 
probability that none of the N particles has entered the 
absorbing sphere. If we assume that the particles are 
randomly distributed in a volume V, then 



Q(t;R) = 



- / drQ 1 {r,t;R) 

V Jr>R 



N 



(2) 



and if we further take the limits N — > oo and V — > oo 
holding the density p — N/V fixed, then this becomes 



Q{t;R) = exp{- P / dr [1 - Q^r, t; R)\ } . (3) 
Jr>R 



The calculation of Qi(r, t; R) for subdiffusive particles 
can be directly adapted from the corresponding calcula- 
tion of this quantity for diffusive particles [l7l fl8L fl9j . 
We introduce the probability density w(r' ,t\r; R) that 
the particle is at location r' at time t if it started at posi- 
tion r at t = 0. As before, an absorbing sphere of radius 
R is centered at the origin r = 0. The survival probabil- 
ity Q\{r, t; R) is related to this probability density by 



Q 1 (r,t;R) = J dr'w(r',t\r;R). 



(4) 



The probability density obeys the fractional diffusion 
equation together with initial and boundary conditions, 



d_ 

at 



w(t' ,t\r; R) = D t 1_7 [DV^w] 



w(r',0\r;R) = <5(r'-r), 
w(R,t\r;R) = 0, 
lim w(r',t\r;R) = 0. 



(5) 

(6) 
(7) 
(8) 



Here V r ' is the Laplacian operator with respect to the 



position r', qD\ 7 is the Riemann-Liouville operator 



D t 7(x,t) = ^7-^-— / dr ■ 



matj {t-ry-r 



(9) 



and I? is a generalized diffusion coefficient. 

Laplace transforming Eqs. (O, 0, and © with re- 
spect to time according to 



g(u)= / dte- ut g(t) 
Jo 

yields for w(r' , u\r; R) 

uw — 5(r' — r) = u 1 'BV^w, 
w(R,u\r-R) = 0, 
lim w(r',u\r;R) = 0. 



(10) 



(11) 
(12) 

(13) 



Note that w(r',u\r; R) is the Green function of Eq. pip . 
The reciprocity of the Green function with respect to the 
arguments r' and r then implies that it also satisfies the 
adjoint equation 

uw - 8{v' - v) = u 1 " 1 DV 2 r w (14) 

along with the boundary conditions [l7| 



w{r',u\R; R) = 0, 
lim w(r' , u\r; R) = 0. 



(15) 
(16) 



Integration over r' then gives for the Laplace transform 
Qi(r, u; R) of the quantity of interest Qx(r, t; R) [IT}, 

uQx - 1 = u^DVlQu (17) 
Q 1 (R,u;R) = 0, (18) 

lim Qi{r,u\R) = -. (19) 

r — >oo U 
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Finally, the inverse Laplace transform yields the evolu- 
tion equation and boundary and initial conditions for the 
survival probability of the particle that started its walk 
at r, 



0_ 



Qi(r,t;R) = [£>V?Qi] 



Qi(R,t;R) = 0, 
Q 1 (r,t = 0;R) = 1, 
lim Qx(r,t;R) = 1. 



(20) 

(21) 
(22) 
(23) 



The last condition is a result of the certain survival at 
any finite time of a particle initially located at r — > 00. 
Due to the spherical symmetry of the problem, the d- 
dimensional Laplacian operator is 



V? = 



1 d 



dr 2 



r dr 



(24) 



and the solution in Laplace space, valid for all d (even 
for non-integer d) is 



uQi(r, u;R) 



-d/2 ifd/2-1 (VrHnjD 



K, 



d/2-1 



^fR 2 ui/D 



(25) 

which will be analyzed in more detail subsequently. Here 
d is the dimensionality and the K's are modified spherical 
Bessel functions of the third kind [20( . Although under 
some circumstances it may be useful and even illumi- 
nating to consider non- integer values of d (for example, 
see |l9j), our results are only physically relevant for di- 
mensions for which the the Laplacian operator (|24p is 
meaningful, that is, for integer dimensions. As an aside, 
we note that this subdiffusive result is related to the sin- 
gle particle survival probability for a normal diffusive par- 
ticle by the relation 

S(t\ 7 )= / drS(T\ 1 = l)T y (r,t) (26) 
Jo 



or, in Laplace space, 



- „7-l 



S(«i f | 7 = l) ) 



(27) 



where S(t\j) is the survival probability associated with 
a physical situation involving a subdiffusive particle and 
Sftl'y = 1) is the survival probability in the same physical 
situation but involving a normally diffusive particle, and 
T 7 (t, t) is the "time-expanding transformation" [2l[ 

T 7 (r,u) = u 7 - 1 exp(-ru 7 ) . (28) 

In our case, Eq. (|2"7|) means that Qi(r,u; R\j) = 
vT'- 1 Qi{r,vri]R\'y = I). [See [U, H2] for more details 
on the integral transformation (|26[) and the scaling rela- 
tion f2?j).] 

Before analyzing these results and using them to obtain 
explicit survival probabilities, we note that one can pro- 
vide an alternative expression to Eq. for the survival 



probability Q(t; R) in terms of Q\. For this purpose, we 
define [3 



f(t;R)= f dr[l-QiM;i?)]- 

Jr>R. 



(29) 



Note that f(0;R) = because Qx(r,0;R) = 1. Taking a 
time derivative of this function, 



(30) 



and using Eq. (|20[1 . we find that 
d 



dt 



R) = - f dr D\-~< [DV 2 r Q 1 (r, t; R)] , 

Jr>R 

= -DvD\-' 1 [ dA-V r Qi(r,t;R), 

J A 



(31) 



or 



r=R 



where Sd = 2n d / 2 /T(d/2) is the surface of a sphere of 
unit radius, and dA is a surface element in the direc- 
tion perpendicular to the surface of the sphere. Since 
/(0; R) = 0, we can write 

f(t; R) = [ dr A/( r; R) = DSdR^Fit; R) (32) 
Jo dr 



with 



F(t; R) = J^dr D^ fcQifa t; R) 



(33) 



and consequently 

Q(t; R) = ex V {-pDS d R d ~ x F{t] R)} . (34) 

We will use this route in our calculations. Our task is 
thus to calculate Qi(r, t; R) and from it F(t; R), to finally 
arrive at the survival probability Q{t; R). 

III. SURVIVAL PROBABILITY: RESULTS 

In this section we present asymptotic results for the 
survival probabilities in all dimensions. 

A. d < 2 

For d < 2 we find that as u — > 

~ T(d/2)(r 2 - d R 2 ~ d ) 

Qi{r,u;R)~ 22 , rfr(2 _ d/2) Ml - 7+7rf/2 > (35) 

and using a Tauberian theorem this leads to 



Qi(r,t; R) 



W 2 ) f ,- d D2 -d, 



2 2 " d r(2-d/2 



\r z - a -R A 



nd/2-1 



(36) 



r(l-7 + 7d/2) 
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But 



B. d = 2 



,1-7 



d 



oA 7 7rQi(r,t;R) 



Or 



r=R 



T(d/2)(2 - d)R 1 - d D d ' 2 - 1 t~ 1 +^ d / 2 

2 2 - d r(2 - d/2)r( 7 d/2) 



(37) 



so that 



F{t-R) = [ dt D]-< ^-Q x ( T \r,R) 
Jo Or 



r=R 



It then follows that 
Q(t;R) ~ exp i -p 



2 d - 1 T{d/2)R 1 ~ d D d / 2 - 1 t^ 2 . 
r(l-d/2)r(l + 7 d/2) 

(47r£>r) d/2 \ 

r(i-d/2)r(i + 7 d/2) I ' 



(38) 



(39) 

Note that this is independent of the radius R of the 
absorbing sphere. Note also that whereas the survival 
probability of the target in the presence of a single 
(sub) diffusive particle decays as a power law with time, 
the decay law becomes a stretched exponential when 
there are many particles at initially random locations. 

With d = 1 this coincides with the result we previously 
reported in (23|, 



-P 



Q(t; R) ~ exp 
When in addition 7 = 1, 

Q(t; R) ~ exp 



r(l+ 7 /2) 



00. 



ipVDt 



(40) 



(41) 



which coincides with the result reported in [19l |. 

Actually, for e? = 1 the survival probability of the tar- 
get in the presence of a single (sub)diffusivc particle can 
be given exactly [and the result (|3"6")) is then the asymp- 
totic expansion of this exact result] as follows: 



1 exp (R - r) y/vTi/D 

Qi(r, u;R) = : 

u u 



so that 

Q 1 (r,t;R) = l-H 1 1 ° 1 



r-R 


(1,7/2)" 




(0,1) 



(42) 



(43) 



Here H is a Fox function. The probability density func- 
tion —dQi/dt of first passage times corresponding to 
Eq. (j42|) a gree s with the one previously obtained by Bal- 
akrishnan [lj] (see also Ref. where —dQi/dt is given 
in terms of a one-sided Levy stable density). When 
7=1 the Fox function reduces to the complementary 
error function and Eq. (|4"3")) then reduces to the classic 
diffusive result 1181 . 



In two dimensions the solution of the fractional diffu- 
sion equation with the appropriate initial and boundary 
conditions leads to 



1 



Qi(r, u;R) — - 



ln{R/r) 



u 7E + ln(i? v /u-r/(4L»)) 
1 \n(R/r) 



(44) 



u \n{Ry/vT</(W)) ' 



0, 



where 7e — 0.577216 is the Euler-Mascheroni con- 
stant. Applying the Tauberian theorem then leads to 
the asymptotic result 

2\n(r/R) _ 2\n(r/R) 

Ql{r ^ R) ~ madFW) = Tw^T- (45) 

To calculate F(t; R) we need to evaluate the fractional 
derivative G(t) = Q D\~ y g{t), with g(t) = l/ln(at). We 
know that 



G(u) = u 1 1 g{u) — lim ■ 



1 



ds- 



(46) 



However, the limit term on the right can easily be shown 



to vanish. Since g(u) ~ [uln(a/u)] 1 as 



0, it 



follows that G(u) ~ [u 7 ln(a/w)] 1 and consequently 
G(t) - t 7_1 /[ r (7) ln(at)] as t -> 00. We thus have that 



F(t;R) = ^- f drG(r) 
Ri Jo 

and consequently 



2V 



i?7 2 r( 7 ) ln(at) 



(47) 



Q(t; R) = exp 



r(l + 7) ln(4£iiT/i? 2 ) 



t — > 00. 



(48) 

While the decay of the survival probability of a single 
particle is thus an inverse logarithm, that of the ensemble 
of particles is a stretched exponential with a logarithmic 
correction. Also, when d = 2 [and also when d > 2 (see 
below)] the survival probabilities do depend on the radius 
R of the absorbing sphere. 



C. d > 2 

For d > 2 we have that 



Qi(r,u; R) 



1 



so that 



i{r,t;R) 



d-2 



u-> 0, 



d-2 



t — » OO. 



(49) 



(50) 
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From this we obtain 
rt 



F(t)=J dt D l -i ^Qi(r,r; 



R) 



r=R 



d-2 



P 



R r(i + 7 )' 

from which it follows that 



S d R d - 2 (d-2) 



Q(t; R) = exp < -p 



r(i + 7) 



oo 



DP 



(51) 
(52) 



t 



(53) 



in agreement with results reported in [l3l [l6j . The sur- 
vival probability of a single particle approaches a con- 
stant at long times, whereas the ensemble survival prob- 
ability decays as a stretched exponential. 

As in the one-dimensional case, for d = 3 it is again 
possible to provide an exact survival probability of the 
target in the presence of a single (sub)diffusive particle 
as follows: 



i(r,u;R) 



1 

u 



so that 

Qi(r,t; R) 



R „W 



Re^~ 




r 


u 


r- R 


(1,7/2)" 




(0,1) 



(54) 



(55) 



This Fox H function reduces to a complementary error 
function when 7 = 1, and in this limit we recover the 
classic normal diffusion result, see e.g. [Hj]. It is inter- 
esting to note that for d = 3 this result is the same as 
Eq. (5.5) in [l| if we take the completely absorbing limit 
fc 7 — > 00 in that formula. Note that the probability den- 
sity function —dQi/dt of first passage times correspond- 
ing to Eqs. (|54|) and (l55f agrees with the one previously 
obtained by Barkai [12 ]. 



IV. CONCLUSIONS 

We have calculated the asymptotic survival probabil- 
ity of an absorbing target of radius R at the origin in 
the presence of one, or of many, (sub)diffusive particles. 
This calculation, which is based on the fractional diffu- 
sion equation, has been carried out for all dimensions, 
with results that in some cases and limits agree with 
known results. Equations (|31)|) and (|3T)|) give these sur- 
vival probabilities for d = 1, Eqs. (j45|) and (|48|) for d = 2, 
and Eqs. (|50|) and ([53")) for d > 3. Thus, in one dimension 
the survival probabilities are respectively of power law 
and stretched exponential form, in two dimensions the 
decay is slower, respectively logarithmic and stretched 
exponential with a logarithmic correction. The result in 
dimensions three [Io| and higher are interesting: the sur- 
vival probability of a single particle goes asymptotically 
to a constant (i.e., the particle may survive forever with 
a finite probability), and this probability is independent 



of the subdiffusive exponent, and thus the same as for a 
normally diffusive particle. The survival probability of 
the target surrounded by a sea of subdiffusive traps does 
decay, again as a stretched exponential that does depend 
on the subdiffusive exponent. Note that in all cases, while 
the mean survival time of the target in the presence of 
a single particle (i.e., the first moment of Qi) is infinite, 
that of a target in a large (infinite) volume containing 
a finite density of particles (i.e., the first moment of Q) 
is finite. This may be an interesting observation beyond 
the biophysical examples mentioned in the introduction, 
for example in the search by an enzyme of a DNA target 
site involving a combination of scanning and relocation 
events in which the relocation times have a power law 
distribution with diverging moments. The mean survival 
time when a single enzyme seeks the target in this three- 
dimensional search diverges (the survival probability of 
the target is a power law without moments), but one of 
many enzymes will reach the target with certainty (the 
survival probability of the target is a stretched exponen- 
tial). 

The quantities that we have calculated lead immedi- 
ately to others frequently used in the literature as well 
as to additional insights. For example, the derivatives 
—dQi/dt and —dQ/dt are the distributions of first pas- 
sage times to the absorbing target. It is also noteworthy 
that the probability that a subdiffusive molecule initially 
at a distance r from the target ever reaches it is the same 
for normally diffusive and subdiffusive particles. In par- 
ticular, since lim^oo Qi(r, t; R) = Q\{r, 00; R) vanishes 
for d = 1 and d — 2, a (sub)diffusive particle will eventu- 
ally reach the target with certainty. On the other hand, 
since for d > 3 this limit is finite, 1 — (R/r) d / 2 , a parti- 
cle escapes the target forever with probability (R/r) d / 2 
whether it is diffusive or subdiffusive. When the tar- 
get is surrounded by a sea of particles, however, since 
lim^oo Q{r, i; R) = Q(r, 00; R) vanishes for all dimen- 
sions, one of the particles, whether diffusive or subdiffu- 
sive, eventually reaches the target with certainty. The ap- 
proach to these asymptotic behaviors of course depends 
on dimensionality and also on the subdiffusive exponent 
7- 

Our immediate future work on this topic will focus 
on the effects of a partially absorbing target, that is, a 
target that does not necessarily "die" upon its first en- 
counter with a (sub)diffusive particle. Some results on 
this case have been reported in the literature, notably 
in [l(| and [ll|, but they do not use the same boundary 
conditions and they also do not consider this situation in 
all dimensions. Clearly, there is still work to be done. 
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